subroutine gcdTest
    use Algorithm_m
    implicit none
    integer, parameter :: a = 60, b = 13
    integer :: d, u, v
    call gcd(a, b, d, u, v)
    write(*, *) d, u, v
end subroutine